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DETERMINATION OF PROPERTIES OF AN OPTICAL DEVICE 



BACKGROUND OF THE INVENTION 

The present invention relates to the determination of properties of an optical 
device under test, e.g. the determination of the group delay of the optical 
5 device. The group delay is a fundamental property of optical devices, such as 
single mode optical fibers or optical components such as Bragg gratings which 
devices are used in the optical transmission of information. 

As for the background of optical transmission of data it has to be said that the 
premier feature of optical fiber is its extremely low loss. This has made it the 

10 dominant transmission medium for long link lengths. The loss characteristics 
of fiber determines were optical communication is practical. At 1550 nanometer 
single mode optical fiber has an attenuation of 0.2 dB/km. This allows fiber 
optic signals to be propagated through very long length of fiber without 
regeneration. Telecommunication systems use the 1300 and 1550 nm windows 

15 for lowest loss in the fiber. Since a telecommunication system must cover a 
very large distance, the aforementioned attenuation of the single strength in the 
fiber is of high importance. Therefore, the loss characteristics of optical fiber 
often limit the distance that a signal can propagate in the fiber. 

However, this is not always the case, in single mode fibers, chromatic 
20 dispersion can limit the distance over which fiber optic signals can propagate. 
Chromatic dispersion describes the fact that the speed of signal propagation 
in the fiber depends on the wavelength of the light. The consequence is that as 
the signal propagates through a long length of fiber, the edges of the waveform 
of the signal start to become more rounded. Eventually, the adjacent bits start 
25 to overlap in time causing the digital waveform to have poor readability. The 
amount of signal rounding depends beyond other parameters on the amount 
of chromatic dispersion in the cable. The problem of this pulse spreading is a 
problem in today's fiber communication because the increasing bit rate of the 
state of the art fiber communication systems, which bit rate reaches numbers 
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of up to 40 Gbit/s per channel, brings the chromatic dispersion in the range of 
the bit resolution of such a 40 Gbit/s communication system. 

To compensate chromatic dispersion in a fiber chirped Bragg gratings were 
developed. In these gratings different wavelengths of the pulse need different 
5 amounts of time traveling through the grating to compensate the dispersion 
caused by the fiber. The concept of this compensation is schematically 
illustrated in Fig. 1. Fig. 1 shows the original pulse 1 that is broadened to a 
spreaded pulse 2 while traveling through a fiber 3. With the help of a circulator 
4 a chirped Bragg grating 5 is introduced in the path of the pulse 2. As 
10 described before the resulting pulse 6 has the original shape of pulse 1 again. 

However, to use such a component as a chirped Bragg grating It is necessary 
to have exact knowledge of the chromatic dispersion caused by such a grating. 
The measurement of chromatic dispersion is accomplished by analyzing the 
group delay through the fiber as a function of wavelength. To measure the 

15 group delay a wavelength tunable optical source such as a tunable laser is 
used to generate coherent light at different wavelengths. In the prior art the 
wavelength of the tunable laser is then incremented step by step and for each 
wavelength step a group delay is determined. Finally, the group delay details 
are used to calculate the chromatic dispersion coefficient The disadvantage 

20 of measuring the group delays step by step is the time, which is needed to 
perform all the wavelength steps and the respective measurements. In other 
words this means it takes quite a long time to get data, which are so precise 
that this precision is sufficient for the needs of the telecommunication industry. 

SUMMARY OF THE INVENTION 

25 Therefore, it is an object of the invention to provide an improved determination 
of a property of an optical device under test, e.g. the group delay of the optical 
device. The object is solved by the independent claims. 

The main advantage of the present invention is that it is possible to measure 
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the group delay and therefore the chromatic dispersion of an optical device, 
such as a fiber or a chirped Bragg grating, in a short amount of time while still 
keeping the precision of the data high enough for the needs of the 
telecommunication industry. This goal is reached by the invention by tuning the 

5 frequency of the coherent light beam of the laser from a maximum to a 
minimum of a given frequency range in a given time interval, e.g. making a 
wavelength sweep through the given wavelength range. Because of the 
inventive correction of effects caused by a non-linearity in the tuning of the 
laser this sweep can be done without detrimental broadening of the resulting 

10 spectrial width. Method and apparatus of the invention avoid the 
aforementioned problems of the prior art and provide for exact data well 
keeping the measurement time low. 

In a preferred embodiment of the invention there is performed a filtering of the 
Fourier transformed first signal with a high pass filter It is further preferred to 
15 use a Manning window as a shape for the high pass filter. Using this filter 
makes sure that a good elimination of not usable data is possible while still 
having important parts of data at the edges of the spectrum within the filter. 
Moreover, it is preferred to use a half of a Manning window as a shape for the 
high-pass filter. 

20 In a further preferred embodiment of the invention this high pass filter is 
adapted to the precision and the shape of the resulting spectrum of the signal 
corrected for a non-linearity of the laser by making an interferometric signal out 
of the corrected first phase signal, Fourier transforming the interferometric 
signal to get a spectral signal, determining a fraction of the maximum of the 

25 spectral signal, determining the abscissas of the intersections of the ordinate 
of the fraction which the curve of the spectral signal, determining the mean 
frequency fmean as the average of the abscissas, band pass filtering the spectral 
signal with a band pass filter having its center at the mean frequency and 
having a width greater than the width of the frequency range. The inventive 

30 concept of adaptive filtering of the signal corrected for a non-linearity of the 
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laser is doubling the inventive success since by eliminating the non-linearity the 
resulting peak of the Fourier transformed spectrum is very sharp so that it can 
be used a filter width which is very much smaller than known filter width in the 
prior art. As can be seen in the following detailed description of the invention 
5 and the respective drawings this inventive concept is resulting in very smooth 
curves showing the group delay. 

In a preferred embodiment the inventive method comprises the following steps: 

* splitting the first initial light beam into a first light beam and a second light 
beam, 

10 • coupling the first light beam into the optical device under test, 

- letting the second light beam travel a different path as the first light beam, 

- superimposing the first and the second light beam to produce interference 
between the first light beam and the second light beam in a resulting first 
superimposed light beam, 

15 - detecting as a first signal the power of the first superimposed light beam as 
a function of time when tuning the frequency of the coherent light beam from 
a minimum to a maximum of a given frequency range in a given time interval, 

- splitting the second initial light beam in a third light beam and a fourth light 
beam, 

20 - superimposing the third light beam and the fourth light beam after each light 
beam has traveled a different path, to produce interference between the third 
and the fourth light beam in a resulting second superimposed light beam, 

- detecting as a second signal the power of the resulting second superimposed 
light beam as a function of time when tuning the frequency of the coherent light 

25 beam from a maximum to a minimum of a given frequency range in a given 
time interval. 
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- using the detected second signal for deriving a non-linearity infomnation about 
a non-linearity in a tuning gradient of the frequency when tuning the frequency 
of the coherent light beam from the maximum to the minimum of the given 
frequency range, and 

5 - using the non-linearity information for correcting effects on the first signal 
caused by the non-linearity to get a corrected first signal. 

It is further preferred to flatten the first and/or the second signal before 
transforming it by approximating it with polynoms of at least second order. 

Additionally it is prefenred to tune the frequency with a tuning velocity of about 
10 5 nm/s. 

Additionally it is prefenred to tune the frequency within a tuning interval of about 
c/1400nm to c/1600nm, c = velocity of light. 

Additionally it is preferred to detect the power at a rate of about 100kHz. 

The term "coherent" in this application means that the coherence length of the 
15 light beam is bigger than the difference of the path of the first and second and 
third and fourth light beams, respectively. 

Other preferred embodiments are shown by the dependent claims. 

It is clear that the invention can be partly or entirely embodied or supported by 
one or more suitable software programs, which can be stored on or otherwise 
20 provided by any kind of data carrier, and which might be executed in or by any 
suitable data processing unit 



BRIEF DESCRIPTION OF THE DRAWINGS 

Other objects and many of the attendant advantages of the present invention 
25 will be readily appreciated and become better understood by reference to the 
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following detailed description when considering in connection with the 
accompanied drawings. The components in the drawings are not necessarily 
to scale, emphasis instead being placed upon cleariy illustrating the principles 
of the present invention. Features that are substantially or functionally equal or 
5 similar will be referred to with the same reference sign(s). 

Fig. 1 shows a chirped Bragg grating in a light path, 

Fig. 2 shows an embodiment of the invention. 

Fig. 3 shows a flow chart of the steps of an embodiment of the inventive 
method, 

10 Fig. 4 shows an illustration of the step 100 of Fig. 3, 
Fig. 5 shows an illustration of the step 106 of Fig.3, 
Fig. 6 shows an illustration of the step 1 16 of Fig.3, 
Fig. 7 shows two possible shapes for filters used in step 1 16 of Fig. 3, 
Fig. 8 shows an illustration of the step 100 of Fig. 3, and 

1 5 Fig. 9 shows the group delay of the grating 1 2 of the example of Fig. 1 . 

DETAILED DESCRIPTION OF THE INVENTION 

Referring now in greater detail to the drawings. Fig. 2 shows a schematic 
illustration of a preferred embodiment of an apparatus 10 for interferometric 
determination of the group delay of a chirped Bragg grating 12 as a device 
20 under test. However, with the present invention it is possible to determine 
optical properties, e.g. the group delay ,of other devices, e.g. fibers or even air, 
also. 

The apparatus 10 comprises a tunable laser 14 as a signal source for a 
coherent laser beam 16, wherein the laser 14 is tunable between 1400 nm and 
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at most 1600nm. The tunable laser 14 can be tuned continuously within a given 
frequency range in a given time interval. The laser beam 16 is coupled into a 
first beam splitter 18, which splits the coherent laser beam 16 In a first Initial 
beam 20 and a second initial beam 22. The first initial beam 20 is coupled into 

5 a second beam splitter 24. With the help of the second beam splitter 24 the first 
initial beam 20 is split into a first beam 26 and a second beam 28. The first 
beam 26 is traveling a different and longer path than the second beam 28, 
since the first beam 26 is coupled via a circulator 52 into the grating 12. 
Reflected by the grating 12 the first beam 26 is coupled via the circulator 52 in 

10 the first path 32 of the first beam 26 again. Then the first beam 26 and the 
second beam 28 are united with a third beam splitter 34 by superimposing the 
first beam 26 and the second beam 28 to a first resulting beam 36. The first 
resulting beam 36 is coupled into a first photodiode 38. The photodiode 38 
transmits its output as a first signal to an analog/digital-converter (ADC) 40 

15 (which can be National Instruments AD-MIO-16DE-10). The ADC 40 is 
connected to a transient recorder 42 that can be a FAST TR1202 of Fast 
Comtec GmbH, Germany. The laser 14 can be a HP81680A of the applicant. 
The transient recorder 42 is coupled to a computer 44 as an evaluation unit of 
apparatus 10. 

20 Second beam splitter 24, first beam 26, second beam 28 and third beam 
splitter 34 built up a first Mach-Zehnder-interferometer. 

Exiting the first beam splitter 18 is a second initial beam 22. A fourth beam 
splitter 46 splits the second initial beam 22 into a third beam 48 traveling a third 
path and a fourth beam 50. The third beam 48 is traveling a different and 

25 longer path than the fourth beam 50, wherein the path length of the third path 
differs from the path length of the fourth path by more than 1 m, preferred by 
less than 20 m, more preferred by less than 9 m, even more prefenred between 
7 and 9 m. This different path is symbolized by some loops 30. The third 48 
and the fourth 50 beam are coupled into a fifth beam splitter 54 which units the 

30 third beam 48 and the fourth beam 50 to a superimposed second resulting 
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beam 56. The secx)nd resulting beam 56 is then coupled into a second 
photodiode 58. The second photodiode 58 transmits its output as a second 
signal to an analog/digital converter ADC 60, similar to the ADC 40. The ADC 
60 is also connected to the transient receiver 42. Fourth beam splitter 46, third 
5 beam 48, fourth beam 50 and fifth beam splitter 54 built up a second Mach- 
Zehnder interferometer. 

The system for determination of the group delay of the optical component 12 
according to the present invention is based on the use of aforementioned 
Mach-Zehnder interferometers, both of which are supplied with light of the 
10 same light source 14. It is necessary to make the determination by 
interferometric measurements since the phase characteristics of the 
component 12 can not be detected directly. 

The calculation of the properties based on the detected first and second signals 
is as follows: 

15 The laser 14 emits light 16 with field strength E dependent on the time t 
according to the following formula: 

Closer (0 =Eo expo- -iTV jy^iOdn (3.1) 

0 

with =-L^ (3.2) 

20 

into the Mach-Zehnder interferometers, Umom denoting the instantaneous 
frequency of the laser 14 . 

In the embodiment of Fig. 2 the frequency of the laser 14 is changed in a linear 
way so that the instantaneous frequency is 

25 v„^(0=v„+?' < (3.3) 

To calculate as a second signal the power of the resulting signal 56 of the 
second Mach-Zehnder interferometer detected with the second detector 58, 
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one puts for the detected voltage Udet that the light of the longer path is delayed 
by T = AL*n/c when arriving at the fifth beam splitter 54. n is the refractive index 
of the medium, AL is the difference of the path length and c is the velocity of 
light. 

5 The detected voltage is determined by the two superimposed fields that travel 
through the different paths of the interferometer: 



10 



I Laser (0+ ^ Laser 

t 

= exp(>.2;r Jv„^(/')iif')+exp(y -In Jv„^(/')^0 

0 0 

t+T t 

l+cos(2^ \v^iV^dr-2n\v^indn 

0 

/+r 



=2 
=2 



(3.4) 



Only regarding of the changing part of this signal and putting r =Ju(t) dt the 
result is: 



U-^ (0=2.cos[2>r.(r„^(/ + r)-r„^(0)]. (3.5) 

Assuming t « t and approximating with differential quotient one reaches: 



15 C/to (0=2-cos 27C ■ -z-j (3. 



6) 



=2 cos[2;r- v^(t) r] (3.7) 



^mom overage over r 



20 The second signal of the second interferometer is therefore ideally, i.e. with a 
linear rise in frequency of the laser 14 over time a cosine, the frequency of 
which is determined by the length AL and the sweep velocity of the wavelength 
over time. 



Because the laser 14 will produce not an ideal sweep over the given frequency 
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range and will produce non linearities in this sweep one can expect a cosine 
which has a frequency which is changing over time. With this second signal it 
is possible according to the present invention to correct the first signal of the 
first interferometer for this non-linearities in therefore getting a resulting 
5 corrected first signal as if the laser 14 had performed a linear sweep. 

In the first interferometer in which the grating 12 is incorporated the group delay 
of this grating 12 is determined. For this purpose as the first signal the detected 
power of the first photodiode 38 is calculated in an analog way as explained 
above. However, in the following the first signal is only determined in its 
10 dependency of the optical frequency u, where as the time dependency is 
determined later on. 

For the detected signal at the photodiode 38 the following formula applies: 
=2[l + cos(>ff^^^(v)Z^^ -fi{v)L)\ (3.8) 

15 

Whereas p is the phase constant of the second path, whereas Pdut is the 
phase constant of the first path through the grating 12. In this formula the fact 
is neglected that the grating 12 can have a reflection coefficient smaller than 
one. However, this has no importance for the changing part of the detected 
20 signal 36. For this changing part it can be written with <|>=P'L: 

t/d;t(v)=2 cos(^^^^(v)-^(v)) (3.9) 
Hilbert-transforming this signal gives the analytic signal: 

i>det(^) = 2-exp(^^^(v)-(^(v)) (3.10) 

25 

and the phase information: 
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10 



<^^(y)=^DVT(y)-'l>(y) (3.11) 

To determine the time dependency of the signal when perfonming a frequency 
sweep of the laser 14 and to find an evaluation mle for the group delay Tg,DLrr. 
the functions <{> and ^dut are now examined. 

For the phase of the path of the beam 28 witiiout the grating 12 which path is 
assumed to be dispersion free it is possible to perform a Taylor approximation 
of the first order: 



^(v) =^(Vo)+-^|v=Vo -(v-Vo) 
a V 

aoi 

=<j>(y^) +r(Vo)-2;r(v-Vo) 

. 2Z,'»(Vo) ^ . 
=(^{y^)+ i-^-2;r(v-Vo) 



(3.12) 



where L is the length of the path without the grating 12. 

The same goes for the path 32 of the beam 26 with the grating 12: 



15 



2 aV ' 

where \|/ is denoting the residual enror. Performing the same steps as in formula 
3.12 one reaches the following formula: 



20 



with 



(3.14) 
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PUT ^ d"^ PUT _ dToyj. 1 ^ dr A" 
dv dX dv dX d_(s^ dX c 

dX\X) 

By putting this formula into 3.1 1 this gives for the phase of the analytic signal: 



^det =^zwr(*'o)-^(n) 

+ 2^[W (Vo)-r(Vo)] (v-Vo) 
dr,, X' (3.16) 
dX c 



=AiP^ (3.17) 

+ 2;iATo (v-Vo) (3.18) 



-^lAo-— ;z:(v-Vo)" (3.19) 
dX c 



-l-T (3.20) 



To determine the group delay of the grating 12 it is sufficient to differentiate the 
20 face <|>detwith respect to oo. This gives: 



dct) 2ji dv In dv dv 

const 

This is the information searched for since only the behavior of the group delay ^ 
25 with respect to the frequency is interesting, rather than an absolute value. 

However, it is easy to extract from formulas 3.17 - 3.20 information about a 
deviation of this linear rise of the group delay with respect to the wavelength for 
gratings with linear chirp. This goes by approximating the measured phase 
signal with a polynomial of second order and subtracting this polynomial of the 
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phase signal. After this only the part 3.20 remains which represents the 
deviation of the group delay with respect a linear chirp grating. 

With the help of the diagram of Fig. 3 the main steps of an embodiment of the 
inventive method are described, whereas some more details of these steps 
5 follow later on: 

In Fig. 3 ovals describe signals and boxes describe measures of the evaluation 
unit 44. First of all the second or reference signal 56 and the first signal 36 are 
towed to the zero line in step 100 (see Fig. 4). The resulting signals 102 and 
104 are now oscillating about the zero line. In step 106 any non-linearities in 
10 the sweep of the laser 14 are corrected with the help of the second signal 56 
(see Fig. 5). Resulting is a signal 108 corrected for any non-linearities in the 
laser sweep. With the help of the aforementioned formulas are evaluated in 
step 110 the group delay 112 and the linear deviation of linear chirp 114. 

With the help of a filter step 1 16 it Is possible to get filtered signals 118 and 120 
15 of the group delay and of the deviation from the linear chirp of the grating 12. 

In Fig. 4 the principle of the function 100 is graphically illustrated. A is the 
amplitude of the signals 36 and 56. 152 denotes the reflection spectrum of the 
grating 12. 124 denotes the ideal center line of the spectrum 122. 

The function 100 is illustrated In Fig. 8 and comprises the following steps: 

20 - making the signal 36, 56 oscillating about a zero line by: 

- determining the points of mean value of the signal 36, 56 by: 

- determining the maximum and the minimum of the signal 36, 56 in a 
predetermined first range of time smaller than the total range of time, 

- determining a mean value between the maximum and the minimum, 

25 - determining the maximum and the minimum of the signal 36, 56 in a 
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predetermined next range of time adjacent the already examined range of time, 

- detemiining a mean value between the maximum and the minimum, 

- repeating the last two steps until the complete time interval is covered, 

- choosing the predetermined range of time by: 

5 - determining the average period of the oscillations of the signal 36, 66, 

- choosing the size of the range so that more than two average periods fit in the 
chosen range of time. 

Figure 5 illustrates a part of the function of step 106 to correct unwitting non- 
linearities in the sweep of the laser 14. The step 106 comprises the following 
10 steps (in the following the symbol T is used instead of " v 

- deriving the non-linearity information by: 

- transforming the second signal 56 to get a Fourier transformed second signal, 

- eliminating the negative parts of the Fourier transformed second signal to get 
a non-negative Fourier transformed second signal, 

15 - retransforming the non-negative Fourier transformed second signal to get an 
analytic signal of the second signal 56, 

- determining the phase of the analytic signal to get as a second phase signal 
the phase as a function of time of the second signal 56, 

- using the second phase signal for determining as the non-linearity information 
20 the frequency as a function of time of the second signal 56, 

- deriving a first phase signal by: 

- transforming the first signal 36 to get a Fourier transformed first signal. 
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- eliminating the negative parts of the Fourier transformed first signal to get a 
non-negative Fourier transformed first signal, 

- retransforming the non-negative Fourier transfomied first signal to get an 
analytic signal of the first signal 36, 

5 - detenmining the phase of the analytic signal to get as a first phase signal as 
a function of time of the first signal, 

- correcting the effects on the first signal 36 caused by the non-linearity by: 

- using the non-linearity information to interpolate the first phase signal of the 
first signal 36 on a linear scale of frequencies to get a corrected first phase 

10 signal, 

For determining the non-linearity infomnation the following steps are perfomned: 

- determining the frequency f(n) of the second signal 56 as a function of n 
discrete points of time, n= 1 N, on the basis of the second phase signal by: 

- detemnining the second phase signal <|>(n) at the n points of time, 
15 - detemnining the maximum <t>max of the second phase signal. 

- using a predetemnined maximum frequency fmax of the frequency range, a 
predetermined average tuning velocity during tuning the frequency and the 
maximum <|)max of the second phase signal to determine for each of the n points 
of time the frequency f(n) according to the fomiula: f(n)=[(fmax-fmin)/»|>nfiax]4>(n). 

20 For the interpolating step, the following suh>-steps are performed: 

- getting the linear scale fiin(n) of frequencies by: 

- using the predetemiined maximum frequency fmax of the frequency range and 
the predetermined minimum frequency fmm of the frequency range to determine 
the linear scale fiin(n) of frequencies according to the fomnula: fiin(n)=[{fmax- 
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fmin)/(N-1)]n, 

- preferably sorting the absolute values of f(n) monotonically, 

- using f(n) for interpolating the first phase signal of the first signal 36 on the 
linear scale of frequencies fiin(n). The Interpolation can be done linearly or with 

5 the help of polynoms or splines. 

The resulting signal 108 which is corrected for any non linearities in the sweep 
of the laser 14 is used in step 1 10 for evaluate the group delay 112 and the 
deviation 114 from the linear chirp of the grating 12. This evaluation is done 
with the help of the above-mentioned mathematical formulas and can 
10 comprises deriving a group delay 1 12 of the optical device under test 12 by 
differentiating the compensated first phase signal with respect to the frequency. 

The method can further comprise the following steps: 

- ignoring at the begin of the tuning a predetermined amount of values of the 
corrected first phase signal to eliminate teething troubles out of the corrected 

15 first signal 108, 

- approximating the group delay 112 with polynoms of at least second order to 
get an approximated group delay, 

- subtracting the approximated group delay from the group delay 1 12 to get a 
non-linear part of the group delay 112. 

20 - using the non-linear parts of the group delay 112 to determine the mean 
signal power of a deviation from a linear of the group delay 1 1 2 of the device 
under test 12. 

- using the square coefficient of the polynomial to determine the mean gradient 
of the group delay 112, 

25 Having evaluated the group delay 112 and the deviation 1 14 it is possible but 
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not necessary to make a filter step 1 16 to filter the signals 112, 114. Actually, 
this filter step is preferably performed before the calculating of the group delay 
112. as follows: 

- adapting the filtering to the shape of the corrected first signal 108 by: 

5 - a: making an interferometric signal out of the corrected first phase signal, 

- b: Fourier transforming the interferometric signal to get a spectral signal, 

- c: determining a fraction, preferred the half, of the maximum of the spectral 
signal, 

- d: determining the abscissas of the intersections of the ordinate of the fraction 
10 with the curve of the spectral signal, 

- e: determining the nriean frequency fmean as the average of the abscissas, 

- f: band-pass filtering the spectral signal with a band-pass filter having its 
center at the mean frequency and having a width greater than the width of the 
frequency range, 

15 - determining the width of the band-pass filter by: 

- a: predetermining or estimating the maximal range GDrange of the group delay. 

- b: determining the mean value GDmean of the group delay according to the 
steps c-e of claim 26, 

- c: calculating the filter width according to the formula: filter width= 

20 fmean(GDrange/GDmean), 

- filtering with a filter having the shape of a Manning window. 
Resulting is a filtered group delay 118 and a filtered deviation 120. 
Figure 6 shows the principle of determination of the center frequency fmean of 
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the interferometric spectmm. 

Figure 7 shows a graph showing different filter shapes for use with the inventive 
nfiethod. a flat top filter shape 130 and a Manning fonn filter shape 140. Figure 
9 shows the group delay 112 over wave length of a chirped Bragg grating 12, 
5 which is a FBGDC-HP1 A of the company JDS Fitel. It is clearly seen that the 
graph of the group delay 1 12 is very smooth with respect to graphs determined 
with methods of the prior art. 
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2. 

15 



20 

3. 



A method of determination of a property of an optical device under test 
(12), comprising the steps of: 

- using a first Initial coherent light beam (20), 

- changing a first initial property of the first Initial light beam (20), 

- coupling the first initial light beam (20) to the device under test (12), 

- detecting a first signal of the first initial light beam (20) received from the 
device under test (12), 

- correcting a non-linearity in the first signal caused by a non-linearity in 
the change of the first initial property by interpolating the first signal on a 
linear scale. 

The method of claim 1 , further comprising the steps of: 

- using a second Initial coherent light beam (22), 

- changing a second initial property of the second initial light beam (22), 

- detecting a second signal of the second initial light beam (22) without 
coupling it to the device under test (12), to discover a non-linearity in the 
second signal caused by a non-linearity in the change of the second initial 
property. 

- using the discovered non-linearity of the detected second signal to 
interpolate the first signal. 

The method of any one of the claims 1-2, further comprising the steps of: 

- producing a coherent light beam (16), 

- splitting the coherent light beam (16) into a first initial light beam (20) and 



-19- 



07-06-2001 



. . , . , |.01 1 13887.2(07-06-2001) 
Agilent T chnologies, InL.. i L 

OurRef.: 20010487 



June 5, 2001 



CLMS 



a second initial light beam (22). 
4. The method of any one of the claims 1 -3. further comprising the step of: 

- detecting the first resulting property simultaneously with the second 
resulting property. 

5 5. The method of any one of the claims 1 -4, 

- changing the first initial property simultaneously with the second initial 
property. 

6. The method of any one of the claims 1-5, 

wherein the first initial property and the second initial property being the 
10 same initial property. 

7. The method of any one of the claims 1-6, 

wherein the initial property being the frequency of the coherent light beam 
(16). 

8. The method of any one of the claims 1-7, further comprising the steps of: 

15 - transforming the first signal in a number of phase signals over a linear 

scale of a number of points of time, 

- transforming the second signal in a number of frequency signals over 
the same linear scale of points of time to discover a non-linearity in the 
second signal caused by a non-linearity in the change of the initial 

20 property, the initial property being the frequency of the coherent light, 

- assigning the transformed first signal to the transformed second signal, 

- interpolating the assigned transformed first signal on a linear scale of 
frequencies. 
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9. The method of claim 8, further comprising the steps of: 

- creating the linear scale of frequencies fiin(n) according to the formula 

fiin(n) = (fmin - fmax) X (n/N), H Element 1 N, N being the number of points 
of time. 

5 10. The method of any one of the claims 1-9, further comprising the steps of: 

- splitting the first initial light beam (20) into a first light beam (26) and a 
second light beam (28), 



- coupling the first light beam (26) to the optical device under test (12), 

- letting the second light beam (28) travel a different path as the first light 
beam (26), 



10 



15 



20 



25 



- superimposing the first (26) and the second (28) light beam to produce 
interference between the first light beam (26) and the second light beam 
(28) in a resulting first superimposed light beam (36), 

- detecting as a first signal the power of the first superimposed light beam 
(36) as a function of time when tuning the frequency of the coherent light 
beam (16) from a minimum to a maximum of a given frequency range in 
a given time interval, 

- splitting the second initial light (22) beam in a third light beam (48) and 
a fourth light beam (50), 

- superimposing the third light beam (48) and the fourth light beam (50) 
after each light beam (48, 50) has traveled a different path, to produce 
interference between the third (48) and the fourth (50) light beam in a 
resulting second superimposed light beam (56), 

- detecting as a second signal the power of the resulting second 
superimposed light beam (56) as a function of time when tuning the 
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frequency of the coherent light beam (16) from a maximum to a minimum 
of a given frequency range in a given time interval, 

- using the detected second signal for deriving a non-linearity information 
about a non-linearity in a tuning gradient of the frequency when tuning the 

5 frequency of the coherent light beam (16) from the maximum to the 

minimum of the given frequency range, and 

- using the non-linearity information for connecting effects on the first signal 
caused by the non-linearity to get a corrected first signal (108). 

1 1 . The method of any one of the claims 1-10, further comprising the steps 
10 of: 

- deriving the non-linearity information by: 

- transforming the second signal to get a Fourier transformed second 
signal, 

- eliminating the negative parts of the Fourier transformed second signal 
15 to get a non-negative Fourier transformed second signal, 

- retransforming the non-negative Fourier transformed second signal to 
get an analytic signal of the second signal, 

- determining the phase of the analytic signal to get as a second phase 
signal the phase as a function of time of the second signal, 

20 - using the second phase signal for determining as the non-linearity 

infonmation the frequency as a function of time of the second signal. 

12. The method of any one of the claims 1-11. further comprising the steps 
of: 

- deriving a first phase signal by: 
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- transforming the first signal to get a Fourier transformed first signal, 

- eliminating the negative parts of the Fourier transiPormed first signal to 
get a non-negative Fourier transformed first signal, 

- retransforming the non-negative Fourier transformed first signal to get 
5 an analytic signal of the first signal, 

- determining the phase of the analytic signal to get as a first phase signal 
as a function of time of the first signal. 

13. The method of any one of the claims 1-12, further comprising the steps 
of: 

10 - correcting the effects on the first signal caused by the non-linearity by: 

- using the non-linearity information to Interpolate the first phase signal of 
the first signal on a linear scale of frequencies to get a corrected first 
phase signal. 

14. The method of any one of the claims 1-13, further comprising the steps 
15 of: 

- determining the frequency f(n) of the second signal as a function of n 
discrete points of time, n= 1, N. on the basis of the second phase 
signal to determine the non-linearity information by: 

- determining the second phase signal ^{n) at the n points of time, 



20 



- determining the maximum <|>max of the second phase signal, 

- using a predetermined maximum frequency fmax of the frequency range, 
a predetermined average tuning velocity during tuning the frequency and 
the maximum (j)max of the second phase signal to determine for each of the 
n points of time the frequency f(n) according to the formula: f(n)=[(fmax- 
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frnin)/<j>max]<t>('^)- 

15. The method of any one of the claims 1-14, further comprising the steps 
of: 



- getting the linear scale fnn(n) of frequencies by: 

5 - using the predetermined maximum frequency fmax of the frequency range 

and the predetermined minimum frequency fmin of the frequency range to 
determine the linear scale fiin(n) of frequencies according to the formula: 

f|ln(n)=[(fmax-fmfn)/(N.1)]n, 

- preferably sorting the absolute values of f(n} monotonically. 

10 16. The method of any one of the claims 1-15. further comprising the steps 
of: 

- using f(n) for interpolating the first phase signal of the first signal on the 
linear scale of frequencies fiin(n). 

17, The method of any one of the claims 1-16, further comprising the steps 
15 of: 

- deriving transmissive and/or reflective properties of the optical device 
under test from the compensated first signal. 

18. The method of any one of the claims 1-17, further comprising at least 
one of the following steps of: 



20 



- deriving a group delay (112) of the optical device under test (12) as a 
function of frequency from the corrected first signal (108), 

- deriving the chromatic dispersion coefficient of the optical device 
under test (12) as a function of frequency from the connected first signal 
(108). 
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19. The method of any one of the claims 1-18, further comprising the steps 

of: 

- deriving a group delay (112) of the optical device under test (12) by 
differentiating the corrected first phase signal with respect to the 
frequency. 

20. The method of any one of the claims 1-1 9, further comprising the steps 
of: 

- ignoring at the begin of the tuning a predetermined amount of values of 
the corrected first phase signal to eliminate teething troubles out of the 
corrected first signal (108). 

21 . The method of any one of the claims 1-20, further comprising the steps 
of: 

- approximating the group delay (112) with polynoms of at least second 
order to get an approximated group delay, 

- subtracting the approximated group delay from the group delay (1 12) to 
get a non-linear part of the group delay (1 1 2). 

22. The method of any one of the claims 1-21 . further comprising the steps 
of: 

- using the non-linear part of the group delay to determine the mean signal 
power of a deviation from a linear group delay of the device under test 
(12). 

23. The method of any one of the claims 1-22, further comprising the steps 
of: 

- using the square coefficient of the polynomial to determine the mean 
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gradient of the group delay (112). 

24. The method of any one of the claims 1-23, further comprising the steps 
of: 

- making the first signal oscillating about a zero line by: 
5 - determining the points of mean value of the first signal, 

- interpolating a curve through these points, 

- subtracting the values of the curve from the first signal to get a corrected 
first signal oscillating about the zero line. 

25. The method of any one of the claims 1-24, further comprising the steps 
10 of: 

- determining the points of mean value by extracting all points with a 
maximum gradient. 

26. The method of any one of the claims 1-25. further comprising the steps 
of: 

15 - making the first signal oscillating about a zero line by: 

- determining the points of mean value of the first signal by: 

- determining the maximum and the minimum of the first signal in a 
predetermined first range of time smaller than the total range of time, 

- determining a mean value between the maximum and the minimum, 

20 - determining the maximum and the minimum of the first signal in a 

predetermined next range of time adjacent the already examined range 
of time, 

- determining a mean value between the maximum and the minimum, 
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- repeating the last two steps until the complete time interval is covered. 

27. The method of any one of the claims 1-26. further comprising the steps 
of: 

- choosing the predetenmined range of time by: 

5 - determining the average period of the oscillations of the first signal, 

- choosing the size of the range so that more than two average periods 
fit in the chosen range of time. 

28. The method of any one of the claims 1-27. further comprising the steps 
of: 

10 - detemnining the points of mean value by: 

- determining the maximum of the Fourier transfomned signal of the first 
signal. 

- using the maximum to detemiine a size of a high-pass filter. 

- filtering the Fourier transformed first signal with the high-pass filter. 

15 29. A method of determination of a property of an optical device under test 
(12). comprising the steps of: 

- detecting a change of a signal with time, being the basis for deriving the 
property, 

- filtering the detected signal by: 

20 - transforming the detected signal to get a Fourier transfonmed signal, 

^ - filtering the Fourier transfonned signal with a filter to get a filtered Fourier 

transformed signal. 
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- retransforming the filtered Fourier transformed signal to get a filtered 
signal, 

- deriving the property on the basis of the filtered signal. 

30. The method of claim 29, further comprising the steps of: 

- the detected signal being corrected for a non-linearity by the method of 
any one of the claims 1 - 28 to get a corrected first signal (108). 

31 . The method of any one of the claims 1-30, further comprising the steps 
of: 

- using the corrected first signal to calculate the corrected first phase 
signal versus frequency, 

- filtering the corrected first phase signal by: 

- Hilbert transforming It before filtering it to get a corrected signal to be 
filtered according to the steps of any one of the claims 29 or 30. 

32. The method of any one of the claims 1-31 , further comprising the steps 
of: 

- filtering the corrected first signal (108) before calculating the group delay 
(112). 

33. The method of any one of the claims 1-32, further comprising the steps 
of: 

- adapting the filtering to the shape of the corrected first signal (108) by: 

- a: making an interferometric signal out of the connected first phase signal, 

- b: Fourier transforming the interferometric signal to get a spectral signal, 

- c: determining a fraction, preferred the half, of the maximum of the 
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spectral signal, 

- d: determining the abscissas of the intersections of the ordinate of the 
fraction with the curve of the spectral signal. 

- e: determining the mean frequency fmean as the average of the 
5 abscissas, 

- f: band-pass filtering the spectral signal with a band-pass filter having its 
center at the mean frequency and having a width greater than the width 
of the frequency range. 

34. The method of any one of the claims 1-33, further comprising the steps 
10 of: 

- determining the width of the band-pass filter by: 

- a: predetemnining or estimating the maximal range GDrange of the group 
delay, 

- b: determining the mean value GDmean of the group delay according to 
1 5 the steps c-e of claim 33, 

- c: calculating the filter width according to the formula: filter width= 

35. The method of any one of the claims 1-34, further comprising the steps 
of: 

20 - subtracting a gradient in the group delay from the group delay. 

- predetermining the maximum range GDrange of the group delay, 

- performing the steps b and c of claim 34. 

- calculating the group delay. 
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- adding the subtracted gradient to the calculated group delay. 

36. A software program or product, preferably stored on a data carrier, for 
executing the method of one of the claims 1 to 35 when run on a data 
processing system such as a computer. 

37. An apparatus of determination of properties of an optical device under 
test, comprising: 

- a first beam splitter (18) in a path of a coherent light beam (16) for 
splitting the coherent light beam (16) into a first initial light beam (20) 
traveling a first initial path and into a second initial light beam (22) 
traveling a second initial path, 

- a second beam splitter (24) in that first initial path for splitting the first 
initial light beam (20) into a first light beam (26) and traveling a first path 
(32) and into a second light beam (28) traveling a second path, 

- a place in that first path (32) for coupling the first light beam (26) to the 
optical device under test (12), 

- a third beam splitter (34) in that first (32) and in that second path for 
superimposing the first (26) and the second (28) light beam after the 
second light beam (28) has traveled a different path as the first light beam 
(26) to produce interference between the first light beam (26) and the 
second light beam (28) in a resulting first superimposed light beam (36) 
traveling a first resulting path, 

- a first power detector (38) for continuously detecting as a first signal the 
power of the first superimposed light beam (36) as a function of time 
when tuning the frequency of the coherent light beam (16) from a 
minimum to a maximum of a given frequency range in a given time 
interval, 
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- a fourth beam splitter (46) for splitting the second initial beam (22) in a 
third light beam (48) traveling a third a path and a fourth light beam (50) 
traveling a fourth path, 

- a fifth beam splitter (54) in that third and in that fourth path for 
5 superimposing the third light beam (48) and the fourth light beam (50) 

after each light beam (48, 50) has traveled a different path, to produce 
interference between the third (48) and the fourth light beam (50) in a 
resulting second superimposed light beam (56) traveling a second a 
resulting path, 

10 - a second power detector (60) for continuously detecting as a second 

signal the power of the resulting second superimposed light beam (56) as 
a function of time when tuning the frequency of the coherent light beam 
(16) from a maximum to a minimum of a given frequency range in a given 
time interval, 

15 - an evaluation unit (44) for deriving optical properties of the optical device 

under test (12), for using the detected second signal for deriving a non- 
linearity information about a non-linearity in a tuning gradient of the 
frequency when tuning the frequency of the coherent light beam (16) from 
the maximum to minimum of the given frequency range, and for using the 

20 non-linearity information for correcting effects on the first signal caused 

by the non-linearity to get a corrected first signal. 

38. The apparatus of claim 37, 

further having computing means capable of performing a least one of the 
further steps according to any of the claims 2-36. 

25 39. The apparatus of anyone of the claims 37 or 38, further comprising : 

a circulator (52) at that place in that first path (32) to enable the apparatus 
(10) to examine reflective optical components also. 
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40. The apparatus of anyone of the claims 37-39, 

where the second beam splitter (24). the third beam splitter (34) and the 
first detector (38) build up a first Mach-Zehnder interferometer. 

and the fourth beam splitter (46) and the fifth beam splitter (54) and the 
5 second detector (58) build up a second Mach-Zehnder interferometer. 
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A method of determination of a property of an optical device under test, 
compnsing the steps of: using a first initial coherent light beam, changing a first 
initial property of the first initial light beam, coupling the first Initial light beam 
5 to the device under test, detecting a first signal of the firet initial light beam 
received from the device under test, correcting any a non-linearity in the first 
signal by interpolating the first signal on a linear scale. 

[Fig. 2 for publication] 
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